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Abstract 

We present a systematic method to derive an ordinary differential equation for any Feynman 
integral, where the differentiation is with respect to an external variable. The resulting differ- 
ential equation is of Fuchsian type. The method can be used within fixed integer space-time 
dimensions as well as within dimensional regularisation. We show that finding the differ- 
ential equation is equivalent to solving a linear system of equations. We observe interesting 
factorisation properties of the D-dimensional Picard-Fuchs operator when D is specialised 
to integer dimensions. 



1 Introduction 



Differential equations are a powerful tool to compute Feynman integrals fTUT^ . Up to now, the 
standard way to derive a set of differential equations for a set of unknown Feynman integrals 
is based on integration-by-parts identities [|8l|9l. This usually results in a system of coupled 
differential equations for a set of so-called master integrals. In addition, there is no guarantee 
that the resulting system of differential equations is of minimal order. In fact, a counter-example 
has been found recently [10]: The two-loop sunrise integral with three unequal masses in two 
space-time dimensions satisfies an ordinary second-order differential equation. Integration-by- 
parts identities lead only to a coupled system of four first-order differential equations. 

We are interested in finding for a given Feynman integral and a given external variable an 
ordinary differential equation of minimal order. In this paper we present an algorithm to solve 
this problem. Our method works for arbitrary space-time dimensions and can thus be used with 
dimensional regularisation. The method reduces the problem of finding an ordinary differential 
equation of order r to the task of solving a system of linear equations. In order to find the 
differential operator of minimal order, one starts with r = 1 . If no solution is found, one increases 
r until a solution is found. The practical limitation of our method is the ability to solve large 
systems of linear equations. 

The method presented here grew out of our previous work on the two-loop sunrise integral 
with unequal masses in two space-time dimensions. In this work we used methods of alge- 
braic geometry and derived an ordinary differential equation for this integral. For the differential 
equation we have to find the differential operator appearing on the left-hand side, as well as the 
inhomogeneous term appearing on the right-hand side. In our previous work this was done in a 
two-step process: We first obtained the differential operator, and then determined the inhomoge- 
neous term by an ansatz. To find the differential operator we used the fact that for this specific 
integral the singularities of the integrand correspond to an elliptic curve. We then computed the 
differential operator from the known cohomology of the elliptic curve. 

For an arbitrary Feynman integral the singularities of the integrand correspond no longer to 
elliptic curves and generalise to algebraic varieties. The cohomology of these algebraic varieties 
is in general not known a priori. The approach of our previous work does therefore not generalise 
trivially to arbitrary Feynman integrals. However, a careful inspection shows that we can encode 
all required unknown information for the differential operator and the inhomogeneous term into 
coefficients of a well-chosen ansatz. Plugging the ansatz into the differential equation will give a 
system of linear equations for the unknown coefficients. We thus reduce the problem of finding 
the differential equation to the (simpler) problem of solving a system of linear equations. 

In addition, our approach has a few bonuses, which we would like to point out: First of all, 
whereas methods of algebraic geometry are often tied to fixed integer dimensions, nothing pre- 
vents us from using an arbitrary space-time dimension within our ansatz. We are thus able to 
derive a differential equation within dimensional regularisation. Secondly, when working within 
dimensional regularisation, we observe interesting factorisation properties of the differential op- 
erator. These factorisation properties are useful when the differential equation is solved term by 
term as a Laurent series in the dimensional regularisation parameter. Thirdly, our ansatz leads 
to a straightforward interpretation of the inhomogeneous term: It can be represented as a sum 
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of Feynman integrals, where one of the original propagators has been removed. These Feynman 
integrals can be considered as simpler Feynman integrals. 

This paper is organised as follows: In the next section we introduce our notation. Section |3] 
presents the method for finding the differential equation. A few examples are discussed in sec- 
tion SI Finally, section [5] contains our conclusions. 



2 Notation and overview 

We start by defining our notation. We consider a Z-loop Feynman graph with m extemal lines 
and n internal lines. We label the extemal momenta by pi, and the / independent loop 

momenta by ki, ki. The momenta flowing through the internal lines are denoted by qi with 
\ < i <n. The momenta qi can be expressed as a linear combination of the extemal momenta pj 
and the independent loop momenta kj with coefficients —1, or 1: 

/ m 

qi = Y.Pijkj+ll^ijPJ^ p,7,a/; e {-1,0,1}. (1) 

The Lorentz invariant quantities in the extemal momenta are 

Sjk = [Pj + Pk] , i<j,k<m. (2) 
The internal masses are denoted by m,- with I < i < n. We call the set 

S = {sjk^mj} (3) 

the set of kinematical invariants. We denote by D G C the dimension of space-time. Within 
dimensional regularisation D is extended from positive integer values to complex values. For 
each internal line the number v,- G N (with I <i<n) gives the power to which the z-th propagator 
occurs. We define the Feynman integral by 
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with V = Vi + ... + v„. The scale ^ is an arbitrary parameter with the dimension of a mass, which 
is introduced to make the integral dimensionless. After Feynman parametrisation we have 



Ig = fo>- (5) 



Eq. ^ motivates the choice of the prefactor in eq. @: The prefactor is chosen such that after 
Feynman parametrisation one obtains a nice formula. The integration in eq. © is over 

A = {[jci :x2 : ... :jc„] G P""^|jc,- > 0,1 < z <«}. (6) 
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(0 is a differential (n — l)-form given by 

n 

CO = ^(-l)-'^^Xj Jjci A... AJxjA... AJ;c„, (7) 

7=1 

where the hat indicates that the corresponding term is omitted. The function / is given by 

/(.„...,.„) ^ (,^)-»/^ (n.;-) ^y^. (8) 

The functions U and are the graph polynomials [fTTI . If one expresses 

j^Xj{-qj + mj) = -££^,M„fc, + £2^,-e,-7, (9) 

j= 1 r=ls=l r= 1 

where M is a Z x / matrix with scalar entries and 2 is a /-vector with four-vectors as entries, one 
obtains 

'li = det(M), T = det{M){-J + QM-'^Q). (10) 

The graph polynomials U and ^ are homogeneous functions in the Feynman parameters xj with 
degrees 

degU = l, degj = / + l. (11) 

The graph polynomial U is independent of the kinematical invariants, while the graph polyno- 
mial ^ depends linearly on the kinematical invariants. Let us denote hy t E S one kinematical 
invariant, t will be our main variable in the following. We seek an ordinary differential equation 
with respect to the variable t for the Feynman integral Iq- We set 

CO, = /CO, (12) 

and hence 



Ig = f^,. (13) 



COf is a (n — l)-form which is homogeneous of degree — (n — 1) in the variables Xj. We look for a 
differential equation of the form 

l('')co, = JP, (14) 

where 



7=0 
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is a Picard-Fuchs operator of order r. The coefficients pj may depend on the kinematical in- 
variants from the set S, the scale /j^, the space-time dimension D and the exponents V,-, but not 
on the Feynman parameters Xi. We normalise the Picard-Fuchs operator such that pr = I. |3 is 
a (n — 2) -form, depending on the Feynman parameters Xi. The differential d is with respect to 
the Feynman parameters jc,. The differential {n — 2) -form (3 is not unique. We may add to P any 
closed {n — 2)-form y without changing eq. (fT4l) . 

Suppose an equation of the form as in eq. (fT4l) exists. Integration yields then 

lW/g = Jd^. (16) 

A 

If the conditions for Stoke's theorem are met, we obtain 

^('•^/g = /p. (17) 
aA 

Within dimensional regularisation there are no obstructions in using Stoke's theorem. If our 
method is used in fixed space-time dimensions, say D = 4, it may happen that although the left- 
hand side is finite, the individual sub-integrals over the faces of A are divergent. In this case one 
replaces P" ^ by a suitable blow-up P and applies Stoke's theorem to the integration in P. 

Eq. (flTl) is the sought-after ordinary differential equation for the Feynman integral Iq. L^''^ is a 
differential operator of order r in the variable t. The right-hand side is given as a sum of integrals 
with (« — 1) Feynman parameters. These integrals correspond to graphs with one propagator less 
and can be considered as simpler. If the "simpler" integrals are known (the right-hand side and 
the boundary values, say at t = 0), then eq. (fTTI) can be used to obtain the Feynman integral Ig- 
This requires to solve the differential equation in eq. (fTTT i. But first, of course, one has to find 
the differential equation. We show that the problem of finding the differential equation can be 
reduced to solving linear equations. 

We comment on our use of the word "Picard-Fuchs". A differential operator of the form as 
in eq. (fT5l) is said to be of Fuchsian type, if all coefficients pj are meromorphic functions of t 
and if pj has at most poles of order (r — j) . Within our method, once the differential equation 
has been found, the Fuchsian property is easily verified a posteriori. The terminus "Picard- 
Fuchs" is usually associated with the following situation: Consider a family of algebraic varieties 
depending on a parameter t and a cohomology class in H"^. One is interested in the variation 
of the cohomology with the parameter t. In order to pass from cohomology classes to functions, 
one considers periods of the cohomology classes. Then the variation of the periods with the 
parameter t is described by a differential equation of Fuchsian type. This equation is refered 
to as a Picard-Fuchs equation. The relation to Feynman integrals is as follows: Fix an even 
integer space-time dimension and suppose that the Feynman integral is finite in this space-time 
dimension. For even integer space-time dimension the integrand in eq. (|5]) is a rational form on 
P"^i with possible singularities on the zero sets of ^ and Zl. In order to resolve the singularities 
one considers a blow-up P of P"^^ and one denotes by Y the strict transform of the variety on 
which the integrand is singular. One further denotes by B the total transform of the integration 
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boundary in eq. (jS]). According to [fT2l the integrand corresponds to a (relative) cohomology 
class of 

H"^{P\Y,B\Br]Y) (18) 

and the Feynman integral to a period of this cohomology class. The variation of this period 
with the parameter t is given by the differential equation in eq. (flTT ). In this way eq. (flTI) is 
a Picard-Fuchs equation. Now eq. (fTTl ) still makes sense if we relax the assumption on even 
integer space-time dimensions and on the finiteness of the Feynman integral. We can consider 
divergent Feynman integrals regulated by dimensional regularisation. Also in this case we find a 
differential equation of the form as in eq. (fTTT ). We continue to call this equation a Picard-Fuchs 
equation. 



3 The method for finding the differential equation 

The method for finding the differential equation can be summarised as follows: In the beginning 
the order r of the differential equation, the coefficients pj (0 < j < r) and the differential form 
P are unknown. We are interested in finding the lowest-order differential equation. To this aim 
we start with r = 1 and an ansatz for the coefficients pj and the differential form |3. From this 
ansatz we obtain a system of linear equations. If a solution for this system exists we are done. 
Otherwise we repeat the game with r + 1 . At least for finite integrals in even integer space-time 
dimensions it is guaranteed that this ansatz works since systems of differential equations coming 
from algebraic geometry are known to have regular singular points. 

Our ansatz is based on the fact, that the singularities of the integrand are given by powers 
of the graph polynomials U and ^. In other words, all singularities lie in Feynman parameter 
space P"^^ on the hyper- surfaces defined by 

11 = 0^ J = 0. (19) 

Polynomials in the Feynman parameters in the numerator may cancel some of the singularities, 
but the main point is that there will be no singularities on new algebraic varieties. This is obvious 
on the left-hand side: Acting with L^''^ on the integrand will only increase the power of ^ in the 
denominator by r. This motivates the following ansatz for the differential {n — 2) -form [3: 

P = i^T'^^^ (20) 

where a is again a (n — 2) -form, but free of singularities. 

Secondly we exploit the fact that Iq is an integral on projective space ¥"^^. This implies 
that the {n — 2)-form a must be homogeneous of degree [(r — 1)(Z + 1) -|- 2] in the variables Xi. 
The degree of homogeneity of a is always a positive integer. From this and the fact that a is 
free of singularities we are led to the ansatz that a can be represented by polynomials of degree 
[(r — 1 ) ( / -f 1 ) -f- 2] in the Feynman parameters xj. 
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Following Griffiths [fT3l we can reduce the degree of the unknown polynomials by one by 
assuming that a can be written in the form 

a = ^ {—\y^^^^\—Xj^aj^+Xj^aj^^^dx\f\...f\dxj^f\.../\dxj^f\...f\dxn, (21) 

where the ai are homogeneous polynomials of degree h= [(r — 1)(Z+ 1) + 1] in the variables Xj. 
For the polynomials a, we assume the most general form. For example, if n = 3 and h = 2 we 
have 



('■) 2 , (0 2 , (0 2 , (0 , (0 , (0 
'^200'''' 1 ^020^2 ' ^002^3 iq-''-1-^2 ~r (^qiiX2X^ -\- G^q^XiX^- 



(22) 



The variables a^'^^ are independent of the Feynman parameters. The most general homogeneous 
polynomial of degree h in n variables has 



n + h-1 
h 



(23) 



monomials. This is the number of possibilities to pick h elements out of a set of n elements by 
not taking the order into account and with repetitions. For a given r we therefore have within 
our ansatz as unknown variables the variables po, pi, Pr-i, which appear in the Picard-Fuchs 
operator (recall that our convention is to take pr = 1), as well as all coefficients appearing in 
the expansion into monomials of the polynomials a,. In the example of eq. (I22l) these are the 
variables a^jl^. In general, we have 

^unknowns ~ r fl 



(r-l)(/+l)+n 
(r-l)(Z+l) + l 



(24) 



unknown variables. The factor n in front of the binomial coefficient comes from the fact that we 
have n polynomials a,. 

Now let us plug our ansatz, specified by eq. (|20l) and eq. (|2T)) into eq. (fT4)) . Working out J(3 
we find 



d^ 



j=l ,-1 J- 



dxi A... A dxj A ... A dxn- (25) 



Plugging this expression into eq. (1141) and comparing the coefficients of dxi A ... Adxj A ...Adx„ 
we obtain for I < j <n the equations 



/ 



i=l 

We multiply both sides in eq. (|26l ) by 

(f?) "^+'^/^"'' ^V-/D/2+r^-v+(/+l)Z)/2+l 



— {-Xiaj+xjai) 



(26) 



n 



-V/+2 



(27) 
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This factor is chosen such that after multiplication both sides of the equation are homogeneous 
polynomials in the variables Xi. We obtain 



n \ I n 



H{ I ri-^^ I {-XiGj + Xjai) + (v,- - 1) I Y[ Xk \ {-Xiaj + Xjai) 



+ (/__v-r+lj [Y[^k\U^^^{-x,aj+Xja,)). (28) 



In eq. (|28] ) we denoted the derivative with respect to t by 

j = ^7. (29) 
dt 

Eq. (|28] ) is our master equation and we should pause a moment to contemplate the main features 
of this equation. First of all, each term of this equation is of degree one or zero in the unknown 
variables (the coefficients pj and the coefficients appearing in the expansion of the polynomials 
a,- into monomials). Secondly, eq. (1281) is homogeneous of degree [n + (Z + 1)(^+ 1)] in the 
variables Xi. Since eq. (|28l) has to hold for all values of the variables Xi, the coefficient c of each 
monomial in the variables xi has to vanish. But each coefficient c of a monomial in the variables 
Xi yields a linear equation c = in the unknown variables. We thus obtain a (possibly large) 
system of linear equations for the unknown variables. In total we obtain by this method 

. - J(^ + l)(^+l)+2n-l\ 

^Vequations - "I I ^^^) 



equations for /Unknowns- The number A^unknowns has been given in eq. (|24|) . Of course, not all 
equations will be independent. With methods from linear algebra we may attempt to solve this 
system. If the system admits a solution, we have found a differential equation for the Feynman 
integral under consideration. In the case where a solution exists, there will be in general more 
than one solution. This is related to the fact that one can always add a closed {n — 2) -form to |3 in 
eq. (fT4l) . This does not affect our method. In order to find the differential equation one can pick 
any solution. 

If the system does not admit a solution, we repeat the exercise by increasing the order r of 
the differential operator L*^''^ by one. In general this will result in a linear system with more 
unknowns and more equations. The practical limitation of this method is the ability to solve 
large systems of linear equations. 

We illustrate our method with a simple example. Consider the integral 

I = ^ j d^xh{\ — x\ — X2) — , 7 = —tXiX2 + {xi+X2)^m^. (31) 

xj>0 
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This corresponds to a one-loop two-point function with equal internal masses in D = 2 space- 
time dimensions. Suppose we search for a first-order Picard-Fuchs operator 



For P our ansatz reads 



= p'^+Po. (32) 
at 



P = y{xia2-X2ai), (33) 



where ai and a2 are homogeneous polynomials of degree 1 in the variables xi and X2. We write 

(1) , (1) 

(2) (2) 

a2 = a\^xi+a}^{X2. (34) 

We have two Feynman parameters and therefore we have to consider the cases j = 1 and 7 = 2 in 
eq. (|28l) . In this example eq. (l28l) is homogeneous of degree 6. The coefficients of the monomials 
in the variables Xi must vanish. This yields the equations 

ni^po + m^ajQ + (2m^ — t) a^^Q — rr^OQ^ = 0, 

2 2 (i) / 2 \ (I) 2 (2) 
m pQ — m a\Q + [Im — t) a^^^ + m a^^ = 0, 

{3n?-t)po + m^a^i^ + 2m^al)\^ + {4m^-t)af(l -m^a'^Q^ = (35) 

This is a system of linear equations in the unknown variables {po,a\Q ,all^^ ,a^^Q ,a^Q^}. In order 
to determine the differential equation we are only interested if the system admits a solution, and 
if this is the case in one particular solution. The system of linear equations for this example is 
easily solved. One finds 



Po 



t{t-Am^) ' 
(1) 2^^np- - 1 [t - Arn^) af^ 



t {t - Am^) 



(2) ^P-rn^ [t - 2m^) + mh [t - Am^) ajj,^ -t{t- 2m^) [t - Arn^) af^ 
""01 = mh{t-Am^) ■ ^^^^ 

The solutions are parametrised by o^q and o'^q . We can choose any values for these two variables. 

In particular we can set a^Q = a^^Q = 0. This freedom of choice corresponds to the fact that we 
can always add a closed form to p. 
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Let us have a closer look at the inhomogeneous term 



(37) 



in eq. (fTTI) . Within dimensional regularisation the integration is over the n faces of the simplex 
A. Let us consider what happens on a particular face. Without loss of generality we can consider 
the n-th face. We consider the restriction of |3 to x„ = 0. If v„ > 1 we have 



Otherwise, if V„ = 1 we find 



PU=o = (-1)"(a'') 



PL„=o = 0. 



(38) 



n 

!=1 



X; 



x„=0 



n-1 



X ^ ( — 1)-' ^ Xj dx\ A ...f\dxj A ... A(ijc„_i. 



Now let 



/ , "mi...m„-^ ■■■-^ 

mj >0,..,m„>0 
mi + ...+m„=(r-l)(Z+l) + l 



(39) 



(40) 



be the expansion of the polynomial a„ into monomials. We recall that a„ is homogeneous of 
degree (r — 1) (Z + 1) + 1 . On the n-th face only the monomials with m„ = are relevant. We thus 
have 



PU=o = (-1)" I 

m\ >0,..,m„_i >0 
mi + ...+ra„_l = (r-l)(/+l) + l 



a-''^ {o(")^ (41) 



with 



An) 

Jm\...m„_\ 



We remind the reader that 



n-l 



, /"-I \ <7;v-(/+l)£>/2 

2NV-/D/2+r-l / -p-,- v,+m,-l \ ^ V ^ W 



(j:y-lD/2+r-\ 



x„=0 



^ (-1)' ^ .x:; <iA:i A ... ^dxj ^ ... ^dXn-\. 



'^lx„=o and J| 



(42) 



(43) 



are the graph polynomials of a graph, obtained from the original one by contracting the edge 
n [fTTll. Therefore 



Jmi...m„-\^ 



(44) 
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p 



Figure 1 : The three-loop banana graph. 

is the integrand of a Feynman integral in (Z) + 2r — 2) space-time dimensions, where the prop- 
agator n has been contracted. Thus the inhomogeneous term in eq. (fTTl ) is given as a linear 
combination of Feynman integrals in (D + 2r — 2) space-time dimensions, where one of the n 
propagators of the original integral has been contracted. We remark that the representation of the 
inhomogeneous term as a linear combination of Feynman integrals is not necessarily unique. As 
previously already mentioned, we can can always add a linear combination of Feynman integrals 
corresponding to a closed form y. 

4 Examples 

In this section we discuss a few examples. Since the size of the system of linear equations 
which needs to be solved grows quite fast with order of the Picard-Fuchs operator, the number 
of loops and the number of propagators, we limit us here to examples which can be solved with 
plain vanilla methods within standard computer algebra systems. A dedicated and more efficient 
method to solve these systems of linear equations is work in progress. 

As already mentioned, the inhomogeneous term is not unique. We can always add a closed 
form to (3, which corresponds to adding a linear combination of (simpler) integrals to the inhomo- 
geneous term, which adds up to zero. As a consequence, the inhomogeneous term usually does 
not have a particular "nice" form and it is usually rather lengthy. On the other hand, the solution 
of the linear system for the Picard-Fuchs operator is unique and usually rather compact. For 
this reason we present here only the Picard-Fuchs operators for several examples. But the reader 
should keep in mind that by solving the system of linear equations one obtains simultaneously 
the solutions for the inhomogeneous terms. 

In particular we would like to give the reader with the chosen examples some indications, 
what order can be expected for the Picard-Fuchs operator. As examples we study the family of 
"banana" graphs. The Z-loop banana graph is a two-point functions with (/ + 1) propagators. As 
an example, the three-loop banana graph is shown in fig. ([T). The /-loop banana graph is defined 
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by 



2W+1-/0/2 



r(/ + l-ZD/2) 



w 



n 



-{p-h 

(/+l)(l-Z)/2) 



l+l 



l+l-lD/2 



(45) 



with 



1 



1i 



and t — p .\x\.D — 2 dimensions these integral simplify to 



xj>0 ^ i=l / 



+ 'U[ Y^Xjtn] 



1 



(46) 



(47) 



and depend only on the graph polynomial ^, but not on the graph polynomial U. In addition we 
consider first the equal mass case mi = m2... = m/_|_i = m. In the equal mass case and in D = 2 
dimensions the Picard-Fuchs operator are particular simple. At one-loop (/ = 1) we find that the 
Picard-Fuchs operator is of order 1 and given by 



(1) 



P 



In? 



d_ 

dt tit-Am^) 



(48) 



This example was discussed in detail in the previous section. For / = 2 the Picard-Fuchs operator 
is of order 2 and reads 



3m^ 



J2 3?2-20m2/' + 9m4 d 
dt^ t{t -m^){t -9nP-) dt t {t - m^) {t - 9ir?) 



(49) 



The results for Z = 1 and 1 = 2 have been known for a long time. With our method we obtain also 
the Picard-Fuchs operator for the three-loop diagram: 



■(3) 



,2^3 



1 



J3 6(^2-1 5mh + 32m^) d^ It^ - 6%mh + 6Am^ d 

1 1 1 

dt^ t{t-Am^){t-l6m^) dt^ {t - Am^) {t - \6m^) dt t^{t-l6m^) 



In D = 2 space-time dimensions the order of the Picard-Fuchs operator for the family of the 
banana graphs is independent of the mass parameters, as long as the masses are non-zero. This 
is related to the fact, that in this case the singularities of the integrand are only determined by the 
graph polynomial ^, but not by U. In the unequal mass case eq. (l48l) generalises to 



P 



d 
dt 



t — m 



t — (mi — mi) 



{nil + mi] 



(50) 
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The generalisation towards arbitrary masses of eq. (@9) has been given in flO]. 

Our method can also be used to derive the differential equation in D dimensions. For the 
one-loop case we obtain in D dimensions 



d ^ {4-D)t^-2{ml + ml)t + {D-2) {m\-mlf 



dt 



2t 



t — {m\—m2) t—{mi+m2 



(51) 



In the two-loop case with unequal masses we find in D dimensions a Picard-Fuchs operator of 
order 4. The increase in the order of the Picard-Fuchs operator compared to the D = 2 case is 
explained by the fact, that away from D = 2 space-time dimensions also the graph polynomial Zl 
can contribute to the singularities of the integrand. However, we observe that this fourth-order 
differential operator factorises in four space-time dimensions into a differential operator of order 
two and two differential operators of order one. We illustrate this first with a simple example: 
Let us consider the integral 



A* 



/ 



d^xd{l —X\ —X2) 



(52) 



This integral corresponds to a one-loop integral in D = 4 space-time dimensions, where one 
propagator occurs to the power two. Unlike I\ in eq. (|50l) . the integral 1\ has a second order 
Picard-Fuchs operator, which can be written in factorised form as 



2^2 I d 



dt 



+ 



'it' -1 {ml + mf) t^ + 5{ml + ml)^t- [mj + mf) [m\-m^^ 



t [t- 



m 



t — (mi — m2j 



(mi +m2) 



d 

dt t 



(53) 

r 



Note that the second factor is independent of the masses. 

Turning back to the two-loop graph with unequal masses in D dimensions we can write the 
fourth-order differential operator L^^^ for this integral in D dimensions as 



(2) 

'reduced 



2\ d (4-D)^(4) 



dt 



(54) 



(2) (4) 

where ^^educed ^'^^ differential operators of order 2 and 4, respectively. Writing out 

these two operators leads to rather long expressions, but the explicit expression is not needed 

(2) 

for the discussion here, ^^educed independent of the dimensional regularisation parameter e = 

(4 — D)/2, while Lg^'* has a Taylor expansion in the parameter e. The factorisation property in 
eq. (|54l ) is very helpful in solving the differential equation: If the differential equation is solved 
term by term as a Laurent series in the dimensional regularisation parameter £, the term with 
the operator has an extra factor of E in front of it and can thus be treated as part of the 
inhomogeneous term. This reduces the problem of solving a fourth-order differential equation 
to solving a second-order differential equations plus two integrations corresponding to the two 
linear factors. 
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5 Conclusions 



In this article we presented a systematic way to obtain an ordinary differential equation of min- 
imal order for a given Feynman integral in a given external variable. The method reduces the 
problem of finding the differential equation to a problem of solving a linear system of equations. 
The method works in D dimensions and can be useful for the calculation of higher-order correc- 
tions in quantum field theory. The inspiration for our method comes from algebraic geometry: 
Feynman integrals can be viewed as periods of mixed Hodge structures, which vary with the 
external variables. This variation is described by a Picard-Fuchs equation. With a suitable ansatz 
we can find the Picard-Fuchs equation, which in turn is the sought-after differential equation for 
the Feynman integral. When working in dimensional regularisation we observe interesting fac- 
torisation properties of the D-dimensional Picard-Fuchs operator when D is specialised to integer 
dimensions. This factorisation property is helpful when the differential equation is solved term 
by term as a Laurent series in the dimensional regularisation parameter. 
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